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We study the renormalization group flow of the average action of the stochastic Navier-Stokes 
equation with power-law forcing. Using Galilean invariance we introduce a non-perturbative approx- 
imation adapted to the zero frequency sector of the theory in the parametric range of the Holder 
exponent 4 — 2e of the forcing where real-space local interactions are relevant. In any spatial di- 
mension d, we observe the convergence of the resulting renormalization group flow to a unique fixed 
point which yields a kinetic energy spectrum scaling in agreement with canonical dimension analysis. 
Kolmogorov's —5/3 law is, thus, recovered for e = 2 as also predicted by perturbative renormal- 
ization. At variance with the perturbative prediction, the —5/3 law emerges in the presence of a 
saturation in the e-dependence of the scaling dimension of the eddy diffusivity. 
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Kolmogorov 's K41 theory [J is the cornerstone of 
current understanding of fully developed turbulence in 
Newtonian fluids. A modern formulation of the theory 
[H is based on the asymptotic solution of the Karman- 
Howarth-Monin equation, expressing energy balance, for 
stochastic incompressible Navier-Stokes equation 
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f-d x P 
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with / Gaussian, incompressible, zero average time- 
decorrelated with correlation 

-< f(x u ti) <8 f{x 2 ,t 2 ) y= S(tn) F(xia) • (2) 
Here -< >- denotes the ensemble average, ® the tensor 



product, Xi 



Xi 



Xj , t ij . t j 



tj and P is a pres- 



sure term enforcing incomprcssibility: d x ■ v = 0. The 
solution of Karman-Howarth-Monin equation predicts in 
any spatial dimension strictly larger than two that the 
energy injected by the external stirring (/) around a 
typical spatial scale L is conserved across an inertial 
range of scales through a constant-flux transfer mech- 
anism, the "energy cascade" , before being dissipated by 
molecular viscosity. In two dimensions, energy and en- 
strophy conservation across the inertial range calls for 
a distinct analysis of the Karman-Howarth-Monin equa- 
tion formalizing the ideas introduced by Kraichnan 
in 0|. The solution predicts a constant flux inverse en- 
ergy cascade from the injection scale towards the fluid 
integral scale. Below the injection scale a constant flux 
enstrophy cascade towards the dissipative scale may take 
place (see e.g. Q). The very existence and properties 
of the enstrophy cascade are, however, sensitive to the 



boundary conditions imposed on ([1} and the eventual 
resence and shape of large scale friction mechanisms 
H [l(| ■ Dimensional considerations based on the solu- 
tion of the Karman-Howarth-Monin equation lead then 
to scaling predictions for statistical indicators of the flow, 
including the —5/3 exponent for the 3d kinetic energy 
spectrum. These predictions convincingly account for a 
wide range of experimental and numerical observations 
(see e.g. 0, [ll[ and references therein). Their first- 
principle derivation is therefore a well-grounded research 
question. A useful tool to pursue this goal is offered 
by the renormalization group, although its application 
to the inquiry of Navier-Stokes turbulence is ridden by 
challenges. Renormalization group analysis [T^-H] can 
be applied only far from the turbulent regime and for a 
very special choice of the random Gaussian field /. This 
latter needs to have in any spatial dimension d a power- 
law spectrum with Holder exponent 4 — 2 e: 



F(x 12 ;m,M) = 

d d p e ip Xl2 T(p) 
(2n) d d- 1 



F(p; m, M) , 



(3a) 
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F(\p; A to, A A/) = A 4 - d - 2£ ^(p; m, M) , (3b) 

with T(p) = I — p ® p/p 2 the transverse projector, 
p :=|| p ||, and m -C M respectively the inverse inte- 
gral and ultra-violet scales of the forcing. The rationale 
for the choice is that for vanishing e the canonical scal- 
ing dimensions of the convective (i.e. dtv, and v ■ d x v) 
and dissipative (i.e. d^v) terms in the Navier-Stokes 
equation tend to the same value. This fact suggests that 
for £ equal zero canonical scaling dimensions may coin- 
cide with the exact scaling dimensions. In this sense, the 
vanishing e case defines a marginal scaling limit around 
which it may be possible to determine scaling dimensions 
by means of a perturbative expansion in e in analogy to 
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what is done for critical phenomena described by a Boltz- 
mann equilibrium (see e.g. [H, EH). For the stochas- 
tic Navier-Stokes equation the situation is, however, not 
conclusive. Renormalization group yields in any spatial 
dimension a kinetic energy spectrum 

£{p) cx P ^-° , 77 2: o = 1 - 4e/3 (4) 

[3 (see also [13 for an exhaustive review) . In the ex- 
ponent labeling emphasizes the possibility of sub-leading 
corrections. Fully developed turbulence in 3d should cor- 
respond to an infra-red dominated spectrum of the stir- 
ring force as it occurs for e > 2. Interestingly, Q recov- 
ers Kolmogorov's result for e equal two. Consistence with 
Kolmogorov theory then requires the exponent in Q to 
freeze for e larger than two to the value —5/3. Within 
the perturbative renormalization group framework, the 
occurrence of such non-analytic behavior can only be ar- 
gued [IH . Direct numerical simulations [l9l . l20j | exhib- 
ited, within a 512 3 -lattice accuracy, a transition in the 
e-dependence of 772 which is consistent with the freezing 
scenario. The situation is, however, completely different 
in two dimensions [U, [22|. On the one hand, pertur- 
bative renormalization group analysis [23j j upholds the 
validity of (|4]) for any e. On the other hand, the asymp- 
totic solution of the Karman-Howarth-Monin equation 
[22l | shows that (0} is always sub-dominant with respect 
to the inverse energy cascade spectrum £(p) oc p~ 5 / 3 , for 
e < 2 i.e. even in the regime where renormalization group 
analysis should apply. Direct numerical simulations up to 
2048 2 resolution give clear evidence of the inverse cascade 
[2ll [22[ . A scenario reconciling these findings may be that 
the Kraichnan-Kolmogorov inverse cascade corresponds 
to a renormalization group non-perturbative fixed point 
which does not bifurcate from the Gaussian fixed point 
at marginality. Evidences of the occurrence of such an 
"exotic" phenomenon, have been given in models of wet- 
ting transitions by "non-perturbative approximations" of 
the Wilsonian renormalization group [24], [25| . More re- 
cently, similar methods gave evidence of the existence 
of a strong coupling fixed point in the Kardar-Parisi- 
Zhang model of interfacial growth [261 ]. yielding scaling 
predictions favorably comparing with direct numerical 
simulations. Motivated by these results, in the present 
contribution we derive the exact renormalization group 
equations for the stochastic Navier-Stokes equation. We 
then investigate them using a "non-perturbative approx- 
imations" similar to the one used in (26|. By this we 
mean, as often done in non-perturbative renormalization 
[27l [28j , truncations of the flow equations based on some 
assumption on the physical properties of the inquired 
system. Specifically, we investigate the consequences of 
the simplest closure compatible with Galilean invariancc 
and with the number of relevant interactions identified 
by perturbative renormalization at small e. The second 
requirement guarantees the existence of a limit where 
the closure becomes exact in the sense that it recovers 
the perturbative renormalization group fixed point. As 
in (2a , |29| , we focus on the exact renormalization group 



equations for the average action or thermodynamic poten- 
tial defined by the stochastic Navier-Stokes equations. 
In striking contrast with the compressible stochastic dy- 
namics studied in [26j, we do not find any evidence of 
a non-perturbative fixed point which may be associated 
to constant flux solutions in general and to the two di- 
mensional inverse cascade in particular. The truncation 
we consider reproduces instead the expected correct scal- 
ing behavior in the regime dominated by real-space lo- 
cal interactions i.e. d — 3 and e < 3/2. Interestingly, 
we observe in any dimension a transition at e = 3/2 in 
the scaling behavior of the eddy-diffusivity. This latter 
deviates from the renormalization group scaling predic- 
tion by freezing from there on in e to its s = 3/2 value. 
This result was previously derived by different methods 
in [301 ] . In spite of the eddy diffusivity saturation, we 
obtain a kinetic energy spectrum scaling in agreement 
to ((H) with no saturation for e > 2. This latter fact is 
not entirely surprising since the particle irreducible ver- 
tices contributing to the approximated renormalization 
group flow, are only a subset of those needed to fully 
reconstruct the flux i.e. the chief statistical indicator 
in Kolmogorov's theory. The structure of the paper is 
as follows. In section [I] we briefly recall the Karman- 
Howarth-Monin equation and its predictions for power 
law forcing. In section [IT] we derive the exact renormal- 
ization group average action for the model. The scope 
of these sections is to provide basic background on tur- 
bulence and functional renormalization to facilitate the 
reading by researchers familiar with one of these subjects 
but not the other. Using the Ward identities imposed by 
Galilean and translational invariance in section 11111 we in- 
troduce our approximations of the exact flow. We write 
the resulting equations in section IIVI where we also out- 
line their qualitative analysis. To simplify the discussion 
we detail auxiliary formulas in appendix [Cj An advan- 
tage of our formalism is that by preserving the structure 
of the exact renormalization group flow it guarantees the 
"realizability" of the "closure" that we impose [3l| . In 
|V]we describe the analytic solution of our equations in a 
simplified limit. Section PVTl reports the result of the nu- 
merical integration of our equations respectively in the 
three and two-dimensional cases. Finally we turn in IVIII 
to discussion and conclusions. 



I. SCALING PREDICTIONS BASED ON THE 
KARMAN-HOWARTH-MONIN EQUATION 

The Karman-Howarth-Monin equation describes the 
energy balance in the putative unique steady-state to 
which Galilean invariant statistical indicators are ex- 
pected to converge. Specifically, if we consider the two- 
point equal time correlation tensor 

C2(a3i2,t)=^ v(xi,t) ® v (x 2 ,t) y , (5) 
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and the three point equal time structure tensor 

S3(asi2) *) =-< Sv (xi2, t) ® Sv (ac 12 ) ® 5« (aji 2 ) >~ , (6) 
<yw(sci2) := «(cci,t) -«(x 2 ,*) , (7) 

a straightforward calculation using incompressibility and 
the inertial range translational and parity invariance 
yields 



d t C+^d x -S-2 K d 2 x C = F , 



(8) 



for C := tr C 2 , F := tr F and S a := S? Ql „ and Einstein 
convention on repeated indices. In any spatial dimension 
strictly larger than two, (|5J) admits an asymptotic solu- 
tion under the hypotheses (see e.g. @ for a detailed dis- 
cussion) that (i-1) statistical indicators attain a unique 
steady state and hence d t C — 0, they are smooth 

for any finite molecular viscosity but (iii-1) the invis- 
cid limit of the energy dissipation exhibits a dissipative 
anomaly 

< -2 lim lim ndlC + -2 Km YmvndlC = .(9) 

kJ.0 ||sb||4.o |H|4,o«4,o 

Under these hypotheses, if the dominant contribution to 
the forcing correlation comes from wave numbers of the 
order m, Kolmogorov's classical result [l| 



lim_ lim 5^ S^ 1 " 2 " 3 (a;) 



1 1 cc 1 1 4.0 k\Q 



2c 



d(d + 2 



-V a {S Pai 6 aaas } , 



(10) 



holds true for m \\ x \\ <C 1, V a being the index cyclical 
permutation operation over a = {ax, a 2 , as). In other 
words, the leading scaling exponent of ^ is 



C3:0 = 1 



(11) 



Dimensional considerations then yield for the kinetic en- 
ergy spectrum scaling exponent the Kolmogorov's scaling 
law 



(12) 



If instead the forcing correlation is a power-law within 
the range of scales M _1 <c|| x ||<C mT l with Holder ex- 
ponent 4 — 2 e, we should distinguish two situations. If 
e < 2, the forcing correlation @ remains well-defined 
in the limit of infinite integral scale mT 1 . In such a 
case (flTjf holds for M _1 <c|| x || <C ^ where we in- 
troduced £ = k /y/F(0), the typical scale below which 
molecular dissipation dominates. Under the present hy- 
potheses £ oc kM 6 ~ 2 , the omitted proportionality fac- 
tor being a dimensional constant independent of n and 
M. This range of scales is not accessible by perturbative 
ultra-violet renormalization group methods. These latter 
may describe instead the range || x || 3> A/ -1 where the 
asymptotic solution of (jHJ states that the leading scaling 
exponent of ^ is 



Dimensional analysis based on (|13|) then recovers the 
renormalization group prediction Q for the kinetic en- 
ergy spectrum. A different scenario occurs for e > 2: 
the forcing correlation has a finite limit if the ultra-violet 
scale M tends to infinity for any finite value of the in- 
verse integral scale m. In the range mT 1 -C || x || <C £, 
I oc nm, £ ~ 2 , (fT0|) holds with possible sub-dominant terms 
with scaling dimension (|13[) . To summarize, the hint 
coming from the Karman-Howarth-Monin equation for 
spatial dimensions d > 2 is that the —5/3 exponent 
stems from the dominance for e > 2 of the constant 
flux over the dimensional scaling asymptotic solution of 
(|8]). Within renormalization group theory two mecha- 
nisms may provide for this phenomenon, either the fixed 
point of the renormalization group undergoes a bifurca- 
tion (see e.g. [32| for an analytically tractable example 
of non-perturbative bifurcation) or irrelevant composite 
operators dominating the limit of vanishing m of the two 
point correlation which become relevant above a thresh- 
old value of e. Examples of such operators are known 
but the scaling dimensions become marginal exactly 
at £ equal two so nothing conclusive can be said about 
freezing. 

The asymptotic analysis of ((8]) in 2d must be treated 
apart in order to take into account enstrophy conserva- 
tion. In particular 4, 5], Kraichnan's theory [7] is epito- 
mized by a more restrictive version of (i-1), which we will 
refer to as (i-2), requiring only Galilean invariant quan- 
tities to reach a steady state. In other words, dtC does 
not vanish. Furthermore, (iii-1) is replaced by a new hy- 
pothesis (Hi- 2) ruling out the occurrence of dissipative 
anomaly for the kinetic energy dissipation: 



lim lim nd 2 C 
f40 ||iB||4.o 



lim lim nd^C 







(14) 



It is worth noting that (iii-2) can be rigorously proved 
to hold true in some setup for the deterministic Navier- 
Stokes [l(| (see also discussion in [HI). We refer 
the reader to 22] for a detailed analysis of the two- 
dimensional Karman-Howarth-Monin equation in the 
power-law case also corroborated by direct numerical 
simulations of ([1]). Here we only summarize the results. 
In the range of scales which can be investigated by pertur- 
bative ultra-violet renormalization group methods, three 
distinct regimes may set in depending upon the value of 
e. For e < 2, the ultra-violet cut-off gives the domi- 
nant contribution to the total energy F(0) oc M 4_2e and 
enstrophy — (d x F)(0) oc M 6 ~ 2e . Correspondingly, the 
inviscid limit in the range M \\ x \\ 3> 1 predicts for the 
leading and sub-leading scaling exponents of ([6]) 



3:0 



1 



& 



3:1 



-3 + 2e 



(15) 



C3:0 



-3 + 2e 



(13) 



This is in agreement with Kraichnan's theory which pre- 
dicts the onset of an inverse energy cascade for wave- 
numbers smaller than the one characteristic of the (to- 
tal) input. The ensuing dimensional prediction for the 
kinetic energy spectrum scaling exponent is (I12[) while 
7/2:1 = 1 — 4/3e only describes a sub- leading correction. 
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For 2 < £ < 3, F(0) oc m 4 ~ 2e and -(fi£f)(0) oc A/ 6 ^ 2£ 
indicate that in the region ?ti~ 1 J> || a; | ^> M -1 the 
third order structure tensor is sustained by an input of 
enstrophy from larger wave-numbers and an input of en- 
ergy from smaller wave-numbers. As a result, the flux 
balances locally in real space with the forcing so that (fT"3|) 
holds true. Finally for e > 3 and in the presence of a 
large-scale hypo-friction [22| both energy F(0) oc m 4 ~ 2e 
and enstrophy— (<9 2 F)(0) oc m 6 ~ 2e inputs are dominated 
by the infra-red mass scale m. As a consequence, a direct 
enstrophy cascade sets in for m || x || <C 1 and 

C 3:0 = 3 & C 3: i = -3 + 2e. (16) 

Again, dimensional analysis based on (| 16[) predicts 

7?2:0 = -3 & Ifrl = 1-4/3 £ , (17) 

with the leading scaling exponent "freezing" at the 
threshold value attained at s = 3. With these results 
in mind, we turn now to the formulation of a non- 
perturbative renormalization group theory with the aim 
of collating scaling predictions for the energy spectrum. 



II. RENORMALIZATION GROUP FLOW FOR 
THE AVERAGE ACTION 

A. Thermodynamic formalism 

For finite infra-red m and ultra-violet M cut-offs of the 
Gaussian forcing ([2]) it is reasonable to assume that the 
generating function 

■ = < e 3M'J+3) y 



J {3,3) 



(18) 



is well defined. The average in (118)) is over the Gaussian 
statistics of the forcing, v(; f+j) is the solution of ([TJ for 
any fixed realization of / shifted by an arbitrary source 
field j, and * denotes the L 2 (M. d x M) scalar product 

3*v{;j):= [ d d xdtj(x,t)-v{x,t;j) . (19) 

JM. d xR 

Functional derivatives at zero external sources of 
(HHJ) yield the expressions of the correlation and response 
(to variations of v with respect to /) tensors of any order. 
The generating function of connected correlations 



(20) 



is equal to minus the free energy of the field theory. In 
particular, with these conventions we have 



/ (3,3) 



$Ja 1 (x 1 it 1 )5j a2 (xi,t 2 ) 



(21) 



3=3=0 



Analogously, the second order response function is 

Sv^jx^h) r w (l,l) 1Q1 ( . \ 



S 2 W(3,3) 



$Jai(xi,tl)5j a2 {x 2 ,t 2 



(22) 



The Legendre transform of the free energy (|20|) specifies 
the average action or the thermodynamic potential of the 
statistical field theory: 

W(u,«) : = sup{j*«+j*«-W( 3)5 )} . (23) 

(3 ,3) 

The Legendre anti-transform of (l23t reconstruct the con- 
vex envelope of the free energy (|2"0)) . In this sense the 
average action may be interpreted as an ultra-violet reg- 
ularization of the theory. The average action is a func- 
tional of the fields (u,u), which are Legendre conju- 
gate to the external sources and which as custom- 
ary will be referred to as "classical fields". As exten- 
sively discussed in [27], [34| the average action provides 
a convenient starting point for non-perturbative renor- 
malization. Dealing with it is conceptually equivalent 
to working with the Wilsonian effective action as done 
by Polchinski in [35j ]. Namely, the corresponding equa- 
tions can in principle be converted in one another by a 
Legendre transform if one identifies the running cut-off. 
The average action offers, as we will see below, some 
technical advantages |28| which significantly simplify the 
formalism. 



B. Flow equations 

A stationary phase approximation to (|18p in the weak 
stirring limit F 4- (see appendix [Bj yields with logarith- 
mic accuracy 

U ~k F ~k u 

Um ~ « * {dt + u ■ da- - k dl)u , (24a) 



d x ■ u — d x ■ u = 



(24b) 



The limit F = describes the trivial steady state of the 
decaying Navier-Stokes equation. We posit that (|24|) pro- 
vides the initial condition for the renormalization group 
flow of the running average action U mr . This flow de- 
scribes the building up of the exact average action U of 
(1181) as a function of an infra-red cut-off suppressing any 
interaction above an infra-red scale m r and recovering 
IA in the limit of vanishing m r . These conditions can be 
matched [13, [36[ if we replace in (TTJ) the molecular viscos- 
ity with an "hyper- viscous" term, local in wave number 
space, 



K I — Y K '.= K -\- K m R ( 



(25) 



with R a function rapidly decaying for large values of 
its argument and diverging at the origin. A convenient 
choice [29| is 



R{p) 



1 



1 



(26) 



3=3=0 



In (|25[) we also introduced the "running" viscosity K mr . 
We will use this extra degree of freedom to constrain the 
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flow to satisfy a renormalization condition on the eddy 
diffusivity. As for the viscosity, we then apply an high- 
pass filter to the Gaussian forcing 



(27) 



such that 

1 

-< /i®/ 2 y=S(t 12 ) X] F W (fi2;mr) , (28) 



i=0 



where we defined 

F( ) (x 12 \m r ) = F (x 12 ; m r , oo) 



(29a) 



trF( ) (p;m r ) = F i 



,4-d-2e 



(rf-l)X(0)(^](?9b) 



and 



F(i) (p; m r ) = F TOr x(i)(p; mr) T(p) , (30a) 

X(i)fe«v) := P e >»r . (30b) 

This latter term describes a local (in the infra-red or for 
m r = O(M)) perturbation of the measure progressively 
suppressed as m r decreases. Locality entitles us to inter- 
pret this term as a renormalization counter-term in the 
sense of [23j, |37J, |38( . Again, we will use the extra freedom 
introduced by F mr to impose a renormalization condition 
on the flow. The replacements (f25j). (|27|) turn (fT8)) into a 
family of generating functions differentiable with respect 
to the parameter m r . A straightforward calculation (see 
appendix I A 1[) yields 

m T .9 mr ^( 3) j) = / d d x\d d x 2 dt 

S 2 Z, 



(m r d mr F) a ^(x 12 ) 



{3,3) 



+ (m r d mr K mr R)(x 12 )dl 



Sj a ^{x ll t)Sj^(x 2 ,t) 

' X2 5j<*i( Xl ,t)5j ai (x 2 ,t) 
Upon defining 

K(x 12 ,t 12 ) := 

Wmr-RCaJia)^! 
_K mr R(a3i 2 )' t F(a?i 2 ) , 



.(31) 



<5(i 



i - / 



(32) 



and 



W\%(x 1 ,x 2 ,t 1 ,t 2 



o (x 1 ,x 2 ,t 1 ,t 2 ) 



(33) 



we can recast pip into the form of an equation for the 
free energy which, in compact form, reads 

itr {(m r d mr K) * (Wgj,, - Wj» , W g>, ,) } . (34) 



Functional derivatives at zero sources of (|34j) spawn a hi- 
erarchy of equations satisfied by the full set of connected 
correlation of the theory. From (|34|) we derive the average 
action flow using the following two observations. First, 
the very definition of Legendre transform (|23p implies 

m r d mr W^^ = -m r d rn M(u,u) ■ 



(35) 

Second, the evaluation of (f33f at zero sources restores 
translational invariance: 

"w( 2 >°) 



W (2) (a; 12 ,i 12 ) 







o (asi2,*i2) • (36) 



The matrix elements of (|36f are specified by the second 
order correlation and response functions ([2"T]) . (|22|) . We 
may refer to them as indicators of the "Gaussian" part 
of the statistics of (fTJ) . We can use (|36l) and the general 
relation 



W (j J) *"(«,«) 



(37) 



following from the Legendre transform (f23|) . to decouple 
the average action into a Gaussian and an interaction 
part [39j : 



I:= W^rW^+M^*] 

(1.1) L {u,u) J 



(j,5) 

Solving this latter relation for W 



(2) 
(3,3) 



(2) 



(2)mti 



/(2) 



(38) 



(39) 



(3,3) 1 

allows us to finally derive the equation for the average 
action : 

1, 



m T d„ 



{il.il) 



- [u , u] * (m r d 7rLr 7Z) -k 



u 


1 


u 





-tr 



(m r d mr TVj 



n=0 



(il.il) ■ 



★ W (2) .(40) 



Some observations are in order. First, the flow equation 
(piO)) is effectively an equation for the reduced average 
action obtained by subtracting the quadratic counter- 
terms associated to the running infra-red cut-off. This is 
desirable because all physical information is indeed con- 
tained in the reduced average action. Second, the flow 
in (|40p does not depend upon the theory under consid- 
eration which instead specify the initial conditions for 
the evolution. This is a formalization of Wilson's idea of 
renormalization group as a flow in the space of the prob- 
ability measures. The fixed point of the flow does not 
depend on the details of the microscopic theory used as 
initial condition for m r = M. It depends instead upon 
the basin of attraction to which the initial condition be- 
longs. Finally, solving (|4U|) exactly is equivalent to solve 
an infinite non-close hierarchy of equations. Perturbative 
renormalization tells us, however, that there are only a 
finite number of relevant coupling, at most two for e <C 1 
and d > 2 [Tt], H3, HI] , determining the scaling properties 
of the stochastic Navier-Stokcs (fTJ . Based on this obser- 
vation, we now turn to the derivation of a truncation 
of the right hand side of (|4"Uf in order to derive explicit 
scaling predictions. 
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III. GALILEAN INVARIANCE AND 
APPROXIMATION 

Perturbative renormalization identifies the number of 
relevant couplings by diagram power counting in unit of 
the ultra-violet cut-off 15]. Relevant couplings corre- 
spond to proper vertices proportional to powers of 
M larger or equal than zero. For the stochastic Navier- 
Stokes equations only IH 1 ' 1 ) for any d and U^ ' 2 -*, for 
d > 2 have non-negative ultra-violet degree. We can 
use this information to hypothesize that pi))) converges 
towards an average action of the form 

W*U (M) *M 



+iu<°' a >(™) 2 + i(u*) 2 U< a ' 1 >*ti 



(41) 



Clearly, the Ansatz closes the hierarchy of equations 
spawned by (|40|) since it is straightforward to verify that 



U 



and by ([23 



(2)int 



(w*U( 2 ^))t 



W (i,i) = u(i,i)t-i 



W (2,0) = _ W (1,D + U (0,2) +W (l,l)t 



Note that 



6v(f + j) 
8 j(x,t) 



(42) 

(43a) 
(43b) 

(44) 



implies that W^ '*) = U (l ' 0) = for any integer i. To 
further evince the rationale behind (1411) we observe that 



Q^fawx \p 2 , u») = (2 n) d+1 

i=l 

2 

x sQ2ui)[tui + Kp(gW(pi,ui)]T(Pi) , (45) 
i=i 

corresponds to a "dressing" of the quadratic coupling in 
(l24l) . Differentiating with respect to p\ at zero wave- 
number and frequency the translational invariant part of 
(|45p provides a convenient non-perturbativc definition of 
the eddy diffusivity. We will therefore refer to (|45|) as the 
"eddy diffusivity" vertex. Also the "interaction" vertex 

U(Pi,wi,p 2 ,w 2 |P3,w 3 ) = 

3 3 



x P (pi,p 2 ) { J (Pi) • T (Ps) ® T(p 2 ) • p 3 } , 



(46) 



(P1.P2) 

admits a similar direct interpretation from (|24p . Finally, 
comparison with (f2~4"f evinces that the "force" vertex 

2 2 

-(2 n) d+1 5& (J2 P t )S(J2 ^) 9 { °' 2) (P 1. wx) T(p x ) (47a) 



,9 ( °' 2) (Pi^i) 



1 1 

^- T ^trF w ( Pl ,TO r )+g(°' 2 )(p 1 , Wl ) , (47b) 

i=0 

describes (minus) the effective forcing correlation. The 
three vertices are, however, not completely independent. 
Galilean invariance constrains the average action to sat- 
isfy the Ward identity (see e.g. poMi^ and appendix IA 2[) 



= r * u 



5 U , a SU 

— * (r • du - r) + — *r ■ du , (48) 

ou ou 



whence it follows after standard manipulations jl5j ] 

O^CPi.wi.O^lp^wa) 

= p 1 d Ul UW(p 1 ,u 1 \p 3 ,u; 3 ) . (49) 



In the context of perturbative renormalization (|49[) is 
used to show that if a parameter fine-tuning ensures that 
(jt 1 . 1 ) i s finite in the limit M tending to infinity so must 
be IK 2,1 ). In general (|49|) is not sufficient to fully specify 
the form of the interaction vertex in terms of iK 1,1 ). If 
we, furthermore, hypothesize 



(2,1) 



1 , 



then ([49]) implies 



g^(puu)=9 {1 ' x) {p) 



(50) 



(51) 



Such an approximation is too rough to give a self- 
consistent model for the full second order statistics. Our 
goal here is more restrictive as it is only to derive self- 
consistent scaling predictions at scales much larger than 
the dissipative. We therefore posit that (HT1) and ([50)1 
may serve for a self-closure able to capture the scaling 
behavior of the zero frequency sector of the theory. We 
also notice that a consequence of imposing (|5U|) . is that 
a generalized Taylor hypothesis Q is verified by the two 
point correlation function for which the dispersion rela- 
tion 



(52) 



holds true. As a final step in the derivation of our ap- 
proximation we rewrite the vertices (|45|) , (|47al) to decou- 
ple explicitly the functional dependence on the cut-off. 
Thus, we couch the eddy-diffusivity vertex into the form 
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(1,1) 



[p;m r ) := 



(1,1) 



(p; TOr) + R[ — 



(53) 



where now -yt 1 ' 1 ) is an unknown non-dimensional function 
which our renormalization group equation will determine. 
Similarly we write 

g^ 2 \p;m r ) := 

[A (0) m 2 - rf - 2£ + A (1) ] p 2 7 (°' 2 )(p;m r ) , (54) 
where we defined the Grashof numbers 

n d f 



A(0) 1 (2 7^4, m 2£ 



(55a) 
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A(i) := 



2-d 



(55b) 



measuring the intensity of the non-local and local compo- 
nents of the stochastic forcing. In the context of pertur- 
bative renormalization the pair (|55l) specifies the running 
coupling constant of the model [13, H3, [HI- In (|55l we 
denoted 



r(f) 



(56) 



IV. APPROXIMATED RENORMALIZATION 
GROUP FLOW 

The Ansatz 

+iu(°' 2 )(*«) 2 + (Tu) • [(Tu) • 0](Tu) , (57) 



with T the transverse projector and (H51 . (|47al) specify- 
ing the Fourier representation of the order two vertices 
summarizes the approximations described in the previous 
section. The insertion of (1571) into the exact renormaliza- 
tion group equation (j40]) yields the equations 



(2)int 



<5u Ql 



j=0 



(1)int 



tr ^ W (2) * 
2 1 <Ht Ql 



-W (2) 



<5W 



(2)int 



(58a) 



j=0 



vd mr {(d-i)5 (0 ' 2) Cp;™r)} 



l(2) 



<5W 



{2)int 



--tr<W"* 

2 



* w 



(2) 



<5W 



(2)tni 



(Sua, 



, f7/(2)mt „ , 

2 fa„, 



(2)int 



Su r 



(58b) 



w=0 



where we defined 

w (2) ;= w (2) + ( TOr .a mr 72.) * W (2) 



(59) 



These equations, the explicit expression of which is given 
in appendix [C] admit a simple diagrammatic interpre- 
tation. Namely if we adopt the symbolic representation 

W (2,0) = _ W (1,1) * y(0,2) ^ [w (l,l)]t = ^ ( g 0a) 

W (l,l) = [lj(M)]t-l _ j ( 6 0b) 



(60c) 



then we can couch equations (I58|) into the form 

m r d mr {/t TOr 7 (1,1) (p; m r )| 



1 



(d-l)p 



■ m 



•9™tr 




, (61a) 



j=0 



m r 9 mr , | [A( ) m 2 d 2 



e + A (1) ] 7 (0 ' 2) (p;^)} 



2 (d - 1) p- 



■m r d m tr 




(61b) 



if we convene to evaluate variations of response (|60a|) and 
correlation (|60b[) lines within the loops according to the 
rules 

?7« «m r -R ( ) + Km r ™> r d mr R ( — ] , (62a) 

\ m r J \ m r J 

m r d mr g( - 2) (p;m r ) « r) F F (i) %i(p, m r ) 
l 

-^f(i) (p • 3 P - d Fi )X(i) (Pi"ir) , (62b) 

where there appear the scaling exponents 

rj K := Tn r —— ln« mr & r/ F := m r — — lnF m (63) 

determined by the fixed point of the renormalization 
group flow and the canonical dimensions 

d Fa =A-d-2e & d Fl = 2 . (64) 



In other words, ([58)1 imply that the functional vector field 
driving the renormalization group flow with our approx- 
imation is obtained by taking the variation of the mode 
coupling equations in a way adapted to (|3"2")l . We sum- 
marize this calculation in appendix [C] Here, we notice 
instead that after turning to non-dimensional variables ( 
P ^ p/ m r) we can rewrite (f58|) as 

[m r d mr - p-d p + ?7 K ]7 (1 ' 1) (p) 

= VF 4 M) (P) - Vn GfrV (p) - (p) , (65a) 

[m r d mr -p - d p + fj F ] 7 (0 ' 2) (p) 

= VF 4°< 2) (p) - Vk (p) ~ Gi°< 2 ) (p) , (65b) 



where 



VF = 



(2 - d- 2e)A( 0) + ??£ A ( i } 
(A( ) + A (1) ) 



(66) 
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The set of the G^'^'s are non-linear convolutions of the 
unknown functions 7' 1 ' 1 ), 7^ ' 2 ) with certain integral ker- 
nels specified by the dynamics. We detail the form of 
these convolutions in appendices IC II and IC 21 In order 
to fully specify the dynamics we need to associate to (|65|) 
two renormalization conditions specifying the coefficients 
(RkJ|) . We require 



A. Fixed point for Am) = Am = 

This is the trivial fixed point. It corresponds to decay- 
ing solutions of the Navier-Stokes equation. 

B. Fixed point for A(o) > 0, Am 7^ 



7 (1 ' 1) ( Po ) = 7 W2J (Po) = l 



(0,2), 



(67) 



where p is the renormalization scale, i.e. the reference 
infra-red scale where we suppose to measure the eddy- 
diffusivity and the force amplitude. Solving the renor- 
malization condition (|67[) for rjp, rj K we obtain 



Vn 



r (lA) M0,2) 


A (1) £,(0,2)" 
A(o)+A(i) *F 


Gi 1 ' 13 


^(1,1) ,-,(0,2) 


A( )+A(i) *F 


[1 + G& x) ] 



Vf 



1 + G 



(1,1) 



r {i,i) r (i,i) ' 

W.F ^*F 



where G\'^' = G^ ,J> {p ) for all k and 



(id) — Q(i,j) ( 
G^ := (p-d p ^)( Po )-G^( Po ), 



q(0,2) 



(2-rf-2e)A (0) 
A(o) + Am 



(68a) 



(68b) 



(69a) 



(69b) 



Ho) -1- A (i) 

The physical motivation behind the renormalization con- 
ditions (157)) is the following. When the running cut-off 
m r is of the order of the ultra-violet cut-off M the aver- 
age action tends to the limit (|23|) with forcing correlation 
dominated by the local component. In such a case we can 
choose 



7 (M) (p)=7 ^ (p) 



(0,2), 



(70) 



for any p: (|70[) indeed specify the initial condition for 
(f6"5| . Irrespectively of m r , we also expect at scales com- 
parable with the integral scale m _1 the bulk statistics 
to be approximately Gaussian, with parameters speci- 
fied by the eddy diffusivity and the renormalized forcing 
amplitude. In between, as m r decreases toward m we 
expect the onset of a non-trivial scaling range in 7 ^ 1,1 - ) , 



(0,2) 



specified by the solution of ([55 
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value of the Grashof numbers A 



(i), 



The initial 
0, 1 parametrize 



the basins of attraction of the truncated renormalization 
group flow. The invariant sets of the planar dynamics 

m r d m ^X( Q ) = -Ami (3r? fc + 2e) 



'•d mr \l) — 



-A(o) (3% • 
-A(i) (3 77 fe 



2 - d - r] F ) 



(71a) 
(71b) 



characterize the possible scaling regimes that our approx- 
imations can capture. A priori we can distinguish four 
cases. 



In such a case the fixed point condition is 
2e 



v. = - T 



rj F = 2-d-2e, (72) 



as predicted by perturbative renormalization [171 l23j . 
Note that negative values of Ai are admissible if the over- 
all "force" vertex remains positive definite. If the corre- 
lation functions also admit a limit as the integral scale 
m tends to zero, we must observe in the scaling range 



and 



2-d-2£ 



(73) 



(74) 



We expect this behavior to be the physically correct for 
< e <C 1 and d > 2. Perturbative renormalization in 
two dimensions 0, |38| also predicts the attainment of 
this fixed point. 



C. Fixed point for A(o) > 0, Am = 
The approximated renormalization group flow equa- 



tions remain well defined in the limit A 



(i) 



0. In such 



a case Gp^\p) — and (|65ap decouples from (|65b[) . 
Furthermore the renormalization conditions yield, self- 
consistently, 



VF =0 . 



(75) 



In other words, the renormalization group equation has 
only one relevant coupling, the eddy diffusivity. This is 
the situation usually faced in perturbative renormaliza- 
tion under the assumption that the spatial dimension is 
bounded away from two. In such a case only U^ 1,1 -* has 
non-negative ultra-violet degree. This implies that there 
is no need to introduce a local counter-term in LK 0,2 ) so 
that F mr is set to zero a priori. The approximated, non- 
perturbative flow here devised reproduces these features. 
It is readily seen that the scaling predictions are then the 
same as in case IIV 51 



D. Fixed point for A 



(o) 



0, A (1) > 



A similar fixed point, if attained, describes an energy 
input dominated by its ultra-violet component indepen- 
dently of e. It is tempting to associate a similar scenario 
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with the 2d inverse cascade. The attainment of such fixed 
point implies 



7] F = 2 - d + 3 rj K 



(76) 



The value of r\ K here needs to be determined dynamically. 

In order to check the realizability of the aforementioned 
scenarios we resorted to the numerical solution of the 
coupled set of equations (|65|) . (|68|) and (|7ip . 



V. A SIMPLIFIED MODEL 

Before turning to the numerical solution of (|65[) . it is 
expedient to analyze a simplified version of the flow. We 
therefore set 



F mr = A(i) — R — 



(77) 



and hypothesize a sharp infra-red cut-off for the power- 
law forcing 



F(p; m r ) = H(p - m r ) F a p 



\-d-2i 



(78) 



Since perturbative ultra-violet renormalization forbids 
non-local counter-terms [37], |38j], these approximations 
are adapted only to the case d > 2. As a consequence, we 
expect (j6"B"j) to converge to the fixed point of section HV CI 



p-a p + |) 7 p)(,) = ^M, (79a) 



(0(0,2), s 

l p .d p -(2-d-e)]^\p) = ^^> 



\(0) 



(79b) 



with Go 1 ' 1 ^, Go°' 2 ' > respectively specified by 
G^\p) C d f'jl-^ 



A ( o) 2p 2 



p2 



[(d-l)p 3 (p-2<f>) + {d-3)p 2 + 2<j>p] 
.g(M)(i)[ g (i,i)(i)+p2 ff (i,i)( P) ] 



(80) 



and 



G { o°- 2 \p) C d f 1 ^{l-f)^ g^{P) 



A ( o) 2 y_r" ^ 4 

[(d - 1) p 2 - 2 dp fc + 2 fc 2 (d + 2 2 - 2)] 



(81) 



5 (M) ( i) 5 (i,i)(p) [ ff (l,l)(l )+ p2 g (l,l) (p)] 
In ([80|) . (|8T|) we used the notation 

P := v / l+p 2 + 2<f>p . (82) 
In the limit p » 1 we can approximate (|79aP as 



Pd p + -r- 7* ' (p) 



(d-1) 



2dp 2 7 . 



(1,1)/ \ (i,i) n x '( 83 ' ) 

* (pn+ (!) 




J - 2 ) 



10 

p 



FIG. 1. (Color online) Result of the numerical integra- 
tion of the dimensionless renormalized functions 7^ 1,1 ^(p) and 
y(°' 2 > (p) for d = 3 and e = 2. The dashed line indicates 
the value 1. Inset: convergence of 7^ 1 ' 1 '(p) from its initial 
(unforced limit) value (p) = 1 toward stationarity (as 

indicated by the arrow). 



whence we infer the leading scaling behavior 
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(1,1) 



(p) 



ptoo J P 3 

under the self-consistence condition 



< e < f 



< £ 



p > 



(d-1) 



2dc +7 P ) (l) (1-f) 



3 

1 1 1. .-: 



> 1 



(84) 



(85) 



Logarithmic corrections may be possible at e = 3/2. Sim- 
ilarly, we can approximate (|79bp as 



(pdp + d ■ 

p 2-d-2e . 



2e- 

(0,2)/ \ 
■7* (P) 



2) 7 i°' 2) b) 



7i M) b) 



P9 p + y 



7i 1 ' 1) b) , (86) 



in the non-dimensional wave-number range defined by 
(|85[). The corresponding scaling prediction is 



.,(0.2)/ 



^ (¥-) i<£ " (87> 

The conclusion is that the model problem kinetic energy 
spectrum should scale in agreement with the prediction 
of the perturbative renormalization group: 

^ 2) (P) 



p 



2-d-2e 



< £ < § 



£(p)~p' 



d _ ip 2-d-2e +7; 



.,(1,1)/ 



p 



(88) 



7. ' ; (p) 

The eddy diffusivity and the force vertices, however, in- 
dividually deviate from the perturbative renormalization 
group prediction. In particular the eddy diffusivity as 
observed first in [30[ saturates to an e independent value 
for £ > 3/2. In Fig [5] we show that the above predic- 
tions compare favorably with the numerical integration 
of (fT9"|) . For < £ < 2, these results are also consistent 
with the direct numerical simulations of [19, [2fJ ■ 



10 



40 



30 



2(1 



10 



FIG. 2. (Color online) Dependence of the fixed point (Ao, Ai) 
(blue dots) on e and d = 3. The fixed point tends toward 
(0, 0) as e 0. 



VI. NUMERICS 

We integrated numerically the set of equations (|65[) 
for the eddy diffusivity and the renormalized forcing am- 
plitude, and equations (|71[) for the coupling constants, 
together with the renormalization conditions 



We proceeded by discretizing the momentum space on 
a logarithmic mesh for p and a linear mesh for 4>. The 
domain of p considered extends from 10~ 4 to 10 4 and 
was covered by 200 points, corresponding to a logarith- 
mic spacing of w 0.092. The angular domain [—1, 1] was 
covered by 100 points. Moreover, in order to improve ac- 
curacy in the scaling range we modeled wave-number log- 
arithmic derivatives p ■ d p using 5-point finite-difference 
expressions. For the mass differential m r d mr we, instead, 
used 2-point finite-difference expressions. This mesh was 
fine enough to observe good continuous convergence of 
the flow equations. 

We integrated the flow equations with initial condi- 
tions (|70p . over m r from 10 down to 10 ~ 9 , using a simple 
Euler explicit method with logarithmic integration steps. 
Finally, we estimated integrals using a trapezoidal rule on 
the linear and the logarithmic mesh. The initial values 
of the Grashof numbers A(o) and Am where randomly 
sampled in the domain 0.01 < A(j) < 10. 

As an example of the numerical integration scheme 
that we have used we show in Fig. Q] the results for 
7 (M) (p) and ^°' 2 \p) for d = 3 and e = 2. Both func- 
tions satisfy smoothly the renormalization condition at 
the infrared limit p = 10 -4 while exhibiting a power-law 
decay in the ultra-violet. We observed the same qualita- 
tive behavior for any values of d = 2, 3 and < e < 4. 
In the inset of Fig. Q] we show the regular convergence of 
the eddy diffusivity toward is final value. 

We first discuss our numerical results in three dimen- 
sions. 
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FIG. 3. (Color online) Scaling exponents A in Eq. [89] of the 
dimensionless renormalized functions 7^ 1 ' 1 '(p) and 7^ 0,2 ^(p) 
and of the energy spectrum £ as a function of e (blue solid 
circles) and d — 3. We also show in red open squares the 
dependence of the scaling exponents 1} K (upper panel) and 
r/F (middle panel) on e. The solid lines correspond to the 
respective renormalization group scaling of equations (|72[) . 
and to p -4e / 3 for the energy spectrum. The dashed lines in 
the upper and middle panels stand for the scalings p" 1 and 
pi-d-4e/3 respectively. 



3d 



For each fixed value of e in (0, 4] , we used the numerical 
scheme described above to integrate the equations 
and obtained, for any initial value of the Grashof num- 
bers for which we found convergence, a single stationary 
solution. This means that for each value of e there exist 
one single fixed point (A(o), \i)) only. In Fig. [21 we show 
the dependence of the fixed point on e. We have noted a 
slower convergence towards the solution as e — > 0, making 
hard to explore the perturbative regime e < 1. Never- 
theless, our results suggest that the trivial fixed point of 
section HV AI is reached in the limit of vanishing e. 

Surprisingly, for small values of e, Am < and be- 
comes positive for a value of e between 1 and 1.25. For 
e > 2, A(o) decreases exponentially, but we always find a 
positive value. 
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FIG. 4. (Color online) Basin of attraction of the fixed point in 
three dimensions, for different values of e. Each yellow (light 
grey) dot, stands for an initial condition of (Am), Am) for 
which convergence was reached. The square with dashed sides 
indicates the domain in which random initial conditions were 
drawn. The light grey solid circles stand for the trajectory of 
the fixed point in the Ami-Am plane and the blue (dark grey) 
circle, to the fixed point for the specific value of e. 



To determine the ultra-violet scaling law as a function 
of e, we computed 



l og7 (^) 



k^^ ee lim 

p->oo lOgp 



(89) 



for (x,y) = (1,1) or (0,2), which defines the scaling ex- 
ponent of the respective function. We denote with A^ £ ) 
the analogous measure for the energy spectrum. 

In Fig. [3] we show the scaling exponents rj K (red open 
squares in the upper panel) and t\f (red open squares 
in the middle panel) as a function of e. Our numerical 
results are in excellent agreement with the theoretical 
predictions (|72[) (solid lines), meaning that our closure 
yields the perturbative renormalization scaling. In the 
same figure we also show the scaling exponent of the di- 
mensionless renormalized functions -yt 1 - 1 ) (upper panel), 
^,(0,2) ^middle panel) and of the energy spectrum (lower 
panel), as a function of e. 

We observe two different regimes. In the first regime, 
for e < 3/2, the eddy diffusivity and the forcing ampli- 
tude scale in agreement with perturbative renormaliza- 
tion, as obtained in (|T3"|) and ([73|l . Instead, for e > 3/2, 
both fields deviate individually from the perturbative 
renormalization prediction. In particular, in this regime 
the eddy diffusivity scales as 7 ( 1 ' 1 ) ~ p^ 1 independe ntly 
of e. This saturation has been predicted first in |30| . 
More interestingly, the deviation of the forcing amplitude 
is such that the energy spectrum scaling is in agreement 
with perturbative renormalization i.e., £ ~ p 1_4e / 3 , for 
all e. Moreover, the deviations of the eddy diffusivity 
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FIG. 5. (Color online) Scaling exponents A in Eq.[89]of the di- 
mensionless renormalized functions 7* ^(p) and 7i°' 2 '(p) and 
of the energy spectrum £ as a function of e (blue solid circles) 
for the d = 3 simplified model of section [V] The red open 
squares in the upper panel correspond to r/ K . The solid and 
dashed lines correspond to the predicted scaling of equations 
((87)) and pg|l. for e < 3/2 and e > 3/2 respectively. 



and the forcing amplitude from the perturbative renor- 
malization coincide with those predicted by our simplified 
model, equations (|84| and (|87|) . 

Finally, we would like to remark some properties of the 
convergence of the numerical scheme that we have used. 
As we mentioned above, the initial seed for the integra- 
tion scheme comprises the initial value of the Grashof 
numbers. We have chosen this initial numbers by drawing 
A(o) and A(i) as random values in the domain [0.01, 10]. 
By doing this, we found that the solution of our numerical 
scheme always converged to the fixed point when e < 3. 
However, for larger e, we noticed that this was no longer 
the case. For e > 3 some of the initial conditions failed 
to converge. This can be seen in Fig. 0] in which we show 
as yellow (light grey) dots, those initial conditions that 
converged to the fixed point. We notice that the basin 
of attraction, limited to the [0.01, 10] x [0.01, 10] domain, 
shrinks as e grows. While we have no ultimate explana- 
tion for this behavior, it may be due either to the very 
small values that A(o) attain for e > 3 or, more trivially, 
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FIG. 6. (Color online) Dependence of the fixed point (Ao, Ai) 
(blue dots) on e and d = 2. The fixed point tends toward 
(0, 0) as e ->• 0. 



to the fact that our numerical scheme fails to converge 
to the fixed point (shown as the blue (dark grey) circle), 
when the initial condition is too far from it. 



B. Single renormalization condition 

We have solved the simplified model of section [V] sim- 
ply by setting r\p = and using the numerical scheme 
described above, by integrating equations (|65]l . (|71a|) and 
(|68al) . In Fig. ([5]) we show the results that corroborate 
the predicted behavior of equations (|84|) . (|87| and (|88|). 

In summary, we have obtained that the stationary so- 
lution to equations (|65*1) is described by equations (|84| . 
(157)) and irrespectively if we impose the system to 
either one or two renormalization conditions. 
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FIG. 7. (Color online) Scaling exponents A in Eq. |89]of the 
dimensionless renormalized functions 7' 1 ' 1 ^(p) and ^°' 2 \p) 
and of the energy spectrum £ as a function of e (blue solid 
circles) and d = 2. We also show in red open squares the 
dependence of the scaling exponents r\ K (upper panel) and tjf 
(middle panel) on e. The solid and dashed lines correspond 
to the predicted scaling of equations (|84p . (|87[) and (|88[) . for 
e < 3/2 and e > 3/2 respectively. 



C. 2d 

In two dimensions the results are in perfect agreement 
with the predictions of equations AMI) . (T57I) and (1551) . 
meaning that the fixed point found is consistent with 
the perturbative renormalization prediction. To start the 
discussion we show in Fig. |6] the fixed point for several 
values of e. The behavior of the fixed point in two di- 
mensions is qualitatively the same as in three dimensions, 
namely the fixed point (A ,Ai) tends to (0,0) as e tends 
to zero; for £ ~ 1, Am < and becomes positive for a 
value of e between 1 and 1.25; for e > 2, A( ) decreases 
exponentially. 

In Fig. [Jj we show the scaling exponent r/ K (red open 
squares in the upper panel) as a function of e, in agree- 
ment with the prediction (fT2"j). Moreover, we also show 
the scaling exponent of the dimensionless renormalized 
functions (upper panel), <y(°< 2 ) (middle panel) and 

of the energy spectrum (lower panel) , exhibiting the same 
behavior as in three dimensions, described by equations 
©, © and (JSFJ. 



Finally, as it was the case in three dimensions, in two 
dimensions we also observed that the basin of attraction 
shrinks for e ^ 3, as is seen in Fig. [5] 

VII. CONCLUSIONS 

The accuracy of the predictions of non-perturbative 
truncations of the renormalization group flow obviously 
depends upon the available physical information about 
the fixed point to which the flow is expected to converge. 
In the case of the stochastic Navier-Stokes equations, we 
used the simplest possible Ansatz for the critical average 
action satisfying Galilean-invariance in order to inves- 
tigate the infra-red behavior of the theory at zero fre- 
quency. This is of course a drastic approximation, cer- 
tainly a priori inadequate to study the multi-time be- 
havior of correlation functions. We were, however, mo- 
tivated by the results of [26[ where it was shown that 
a similar approximation appears to be able to capture 
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calness" of interactions. This means interactions which 
are local not in real space but in wave-number space. 
Localness of interactions was the motivation behind the 
attempts to derive scaling in turbulence models by imple- 
menting renormalization group transformations through 
progressive resummation of the infra-red degrees of free- 
dom (see e.g. [H, 13]). Our results indicate that further 
explorations of non-perturbative renormalization of the 
stochastic Navier-Stokes equation in the localness of in- 
teraction regime should necessarily encompass the full set 
of proper vertices contributing to the flux. As a conclu- 
sive remark we observe that non-perturbative renormal- 
ization methods may have spin-offs for engineering ap- 
plications: in three dimension for < e < 2 the results 
of the present paper already provide a simple model for 
the computation of eddy-diffusivities which can be used 
in large eddy simulations (45| . 



FIG. 8. (Color online) Basin of attraction of the fixed point in 
two dimensions, for different values of e. Each yellow (light 
grey) dot, stands for an initial condition of (A(o),A(i)) for 
which convergence was reached. The square with dashed sides 
indicates the domain in which random initial conditions were 
drawn. The light grey solid circles stand for the trajectory of 
the fixed point in the A( )-A(i) plane and the blue (dark grey) 
circle, to the fixed point for the specific value of e. 



the existence of non-perturbative fixed point for the KPZ 
stochastic partial differential equation. This latter model 
shares with the stochastic Navier-Stokes equation invari- 
ance under Galilean transformations and convergence to- 
wards a non-Boltzmann steady state. An important dif- 
ference between these two models resides, however, in 
the non-locality of the interactions that the incompress- 
ibility condition brings forth. In retrospective, this dif- 
ference explains why our approximation can reproduce 
and extend the predictions of the perturbative renormal- 
ization group to a finite range of values of the Holder 
exponent. On the one hand, direct numerical simula- 
tions corroborate these predictions in the physical range 
(0 < e < 2 for d > 2 and 2 < e < 3 in d = 2) 
where phenomenological insights support the idea that 
the critical theory should be described by an average ac- 
tion with relevant local (in real space) interactions. On 
the other hand, in the regime of fully developed tur- 
bulence, all evidences point at the fact that the local- 
interaction energy spectrum (|4]) becomes a sub-leading 
correction to the constant- flux —5/3 kinetic energy spec- 
trum. Our approximations fall short of explaining this 
behavior. There are two types of mechanisms which may 
account with renormalization theory for the constant flux 
scaling predictions. Either the dominance of statistical 
indicators by composite operators perturbing the scaling 
relations dictated by the local average action fixed point 
(see discussion in [17|) or the onset of a non-perturbative 
fixed point, the scenario that we tried to explore in the 
present contribution. Fully developed turbulence is be- 
lieved (see discussion in Q) to be characterized by " lo- 
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Appendix A: Variations of the generating of function 

1. Renormalization group flow 

Let us consider the deformation of (fTJ) induced by the 
replacements k h- K + K„ lr R and / h-> f'+j. We suppose 
that / is obtained from applying to / an high pass filter 
with infra-red cut off m r . We have then 

■m r d rrir Z {3 - 3) =~< e 3 * v j* (m r d mr v) >- , (Al) 

with 

m r d m7 .v{x, i; j + f) = — * 

{(m r d mr K mr R)-kd 2 v + (m r d mr f)} . (A2) 
In (|A2|) the fluctuating response function satisfies 
<fo(jci,fi) 



Sj(x 2 ,t 2 



Vi 2 <ii 



(A3) 



We furthermore interpret the product of the time 5- 
correlated Gaussian field / with other functionals in 
(|A1[) according to Stratonovich convention in order to 
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preserve ordinary calculus. Using (|A3[) we can write 
-< e 3 



$3 

5 e 3 * v 

=~< — * { m rd m K m R) * d 2 v >~ 

= tr(m r d mr n mr R) * . (A4) 

Furthermore, a functional integration by parts yields 



-< e 



* {m r d mr f) >- 



1 6 2 e 3 * v 

= - -< (m r d mr F') * — !>- , (A5) 

2 ojdj 

the factor 1/2 being a consequence of Stratonovich con- 
vention. 



2. Ward identity 

Let r t : K ->■ M d a smooth path. The generalized 
Galilean transformation 



a; = a; + e r t , 



(A6a) 



v = v + eft , (A6b) 

leaves (JTJ) invariant in form when if accompanied by the 
redefinition of the forcing / = / + er*. We must have 
therefore 



z { 3 r 3 ) ~ Z b,o) 



(A7) 



When we differentiate this equality at e equal zero and 
use (|A2|) we obtain after standard manipulations (see e.g. 



= f * 



6W, 



{3,3) 



j* [r ■ d- 



SW, 



(3,3) 



+ 



S3 



f J + J* [r ■ d 



SW, 



(3,3) 



$3 



• (A8) 



An alternative way to derive the results of this appendix 
is based on the Janssen-De Dominicis Eg, path inte- 
gral representation of (|18|) . We refer to [42J for a detailed 
presentation. 



Appendix B: Janssen— De Dominicis path integral 
and optimal fluctuation 

The Janssen-De Dominicis [HI, [47j representation is 
the formal measure on path space obtained by requir- 
ing through an infinite dimensional product of Dirac 5- 
functions that at any space-time point ([TJ be satisfied. 
The resulting expression is then averaged over the real- 
izations of the stochastic forcing. We obtain 

Z(3,3) = J D[v]D[v]e~ A , (Bla) 



A = 



v * F * v 



-iv -k \{p t - nd x )v + J(v ■ d x v) - 3} 



(Bib) 



A precise meaning to (|B1I) can be given on a space- 
time lattice using a pre-point discretization dt (v ■ dtv) ~ 
v(U) ■ [v(t l+1 ) - v(ti)}, dt /(©(<), »(<)) ~ dt f(v(U), v(U)) 
for all other terms in (IBlbl) . Notice that in the limit of 
vanishing stirring F 1. 0, (|B1I) recovers the Fourier repre- 
sentation of a product of Dirac ^-functions localizing the 
measure over the deterministic decaying dynamics. In 
this sense (|B1[) remains meaningful also as a formal mea- 
sure inclusive of compressible fluctuations. From (|Blb|) 
a stationary phase approximation yields the weak noise 
limit of the free energy W( Ji j) around an optimal fluctu- 
ation v* . As usual [48}, the stationary phase condition is 
derived by closing a contour in the complex variables 



V = Vsr + IV cj 



(B2) 



which decomposes (IBlbl) into the real and imaginary 
parts 

vx + F + vm 
^ A (3-j) = ^ 3 + 

«s* j(<9 t - ndl)v + T(v ■ d x v) - ^F*v a - j j ,(B3a) 
^■^■(3,3) = 

-vsr *{(d t - k&1)v + T(v ■ x v) - F * v Q - 3} .(B3b) 

The stationary phase condition ^sA(j .j) = can then be 
solved for vq and leaves with a convex functional of the 
principal field v. Assuming that we can minimize such 
functional for some assigned boundary condition, we find 
within logarithmic accuracy 



|| (d t -Kd x )v*+T(vd x vy-3\\l- 



(B4) 



where || v ||p stands for || v ||p= v*F 1 *v. The Legendre 
transform gives the conditions 



u 



(B5a) 



u = F" 1 * {{dt - Kdl)v* + T(v ■ d x v)* - 3} (B5b) 
whence we finally obtain 

U ~ u-k (d t + u ■ d x - k d x )u — i« * F * u . (B6) 

Appendix C: Explicit expression of the convolutions 

An alternative derivation of the renormalization group 
equations is obtained if we observe that we may interpret 
the free energy defined by the Ansatz for the average 
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action ([57)1 as solution of a formal Janssen-De Dominicis 
[M S3 path integral 



W(j,j) = lime In / D[u]D[u]e 



Computing the right hand side in a perturbative expan- 
sion in powers of the interaction vertex (j46 j) . (|50[) we ob- 
tain by standard diagrammatic techniques 



«m r P 2 7 (1,1 HP/ ?7 v) = 

d d k (1-0 2 ) N^(p,k,q>)g^ 2 Hk) 
(2ir) d 2g(^)(k)D 1 (p,k,<f>) 



(C2) 



and 



[A(o)^- 2E + A (1) bV ' 2) (]'K) = 
d d k (l - </> 2 ) A^ (0 ' 2) (p, k, <j>) g(°- 2 > (g) g (°- 2 ) (fc) 
(2^F 4 5 (i.D(fc) 3 (i.i)(Q)^ 1 (p,fc, ( /)) 



(C3) 



We recover equations (|65j) by taking the logarithmic 
derivative m r d mr of both sides of (|C2[) . (IC3[) Note that 
in ([Cl]). El we denoted 



Q := p k 



(C4) 



and </> the cosine between the external p and the integra- 
tion k wave-numbers: 



p ■ k 

p k 



(C5) 



We also defined the auxiliary integrand factors 

D 1 (p,k^)=k 2 g^ 1 \k) + Q 2 g^ 1 \Q) , (C6) 

C 2 (p,fc,0) = 2fc 2 .g( 1 ^ 1 )(fc) + Q 2 3 ( 1 ' 1 )(g) , (C7) 
and the constants 



Gt =(<*"!) / ^(W 2 ) 



(C8) 



Finally, the convolutions depends upon certain integral 
kernels which stem from the expansion up to one loop 
accuracy of the Ansatz average action (|57p . These are 



(d - 1) p 3 {p - 2 k) + k 2 p [(d - 3) p + 2 <f> k] 
k 2 (p 2 + k 2 -2kp<j)) 

N (1A Hp,k,cb) := 
pk [{d - l)pk - 2 (p 2 + k 2 - 2pk<f) <f>] 
k 2 (p 2 + k 2 -2kp<t>) ' 



, (C9a) 



(C9b) 



for the eddy diffusivity vertex, (|C9b[) will be needed be- 
low, and 



JV<°' 2 >(p,M) := 

p 2 [(d- l)p 2 - 2dpk<j) + 2k 2 (d + 2(j) 2 - 2) 
k 2 (p 2 + k 2 - 2kp4>) 2 



(CIO) 



for the force vertex. Finally in (|65[) there appear terms 
of the form 



.(CI) G?' j) (p) 



?± r * fc d / # a - ^ vpn (p, fc, ^) , (on) 



2p 2 7 Q fc 



with Z taking values {F, «, o} and Vg (p, fe, </>) the non- 
linear convolutions specified below. 



1. Equation for the eddy diffusivity vertex 

The following three non-linear convolutions enter 
(I65a| : 



v£ u) (p,M) := 
with coefficient 77^ 

v; (m) (p,m 



iy( 1 ' 1 )(p,fc, ( />)A (1)X (i)(fc) 
g(^)(k)D 1 (p,k,4>) ' 

_ j?(fc) 
:_ [^(p.M)] 2 * 
D 2 (p, fc, </>) N^'^ip, k, 4>) g^ 2) {k) 

fc 4 A>(M)(p ; fc^) g (0.2) ( Q) | 

QV^KQ) J ' 

with coefficient ry K , and 

y o ( 1 - 1 )( P ,fc,0) = 

iV^ 1 ) (p, fc, 0) Ei=o A M (* • flfc - <*f w )xw (*) 



(C12) 



(C13) 



g^{k)D x {jp,kA) 

D 2 (p,k,cl > )N^(p,k^)g^ 2 Hk) 



(k-d k R)(k) 
[A(p,fc,0)] 2 

k i N^' 1 \p,k,4>)g ( - ^(Q) 



(C14) 



QVi,D(Q) 
with coefficient equal to the unity. 

2. Equation for the force vertex 

The following three non-linear convolutions enter 
((65b]) : 

^ F (Q) (A;) -Di (p, fc, 0) ' (C15) 

with coefficient t?f, 

Fi°' 2 )(p,fc,0) :=iV^ 2 )(p,fc,^)x 
g(°- 2 )(Q)g^ 2 \k)R(k)D 2 (p,k, 



g(^)(Q)[g(^)(k)] 2 [D 1 (p,k,< 



(C16) 
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with coefficient r/ K , and with coefficient equal to the unity. 

\ H A « { k ' d k- d F{i) ) X(i) (k, n) 
(k-d k R)(k)g^(k) D 2 (p,k,4 



gW(k) I>i(p,M) 



(C17) 
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